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Abstract 

The nonlinear evolution of a unstable electrostatic wave is considered for a 
multi-species Vlasov plasma. From the singularity structure of the associated 
amplitude expansions, the asymptotic features of the electric field and distri- 
bution functions are determined in the limit of weak instability, i.e. 7 — > + 
where 7 is the linear growth rate. The asymptotic electric field is monochro- 
matic at the wavelength of the linear mode with a nonlinear time dependence. 
The structure of the distibutions outside the resonant region is given by the 
linear eigenfunction but in the resonant region the distribution is nonlinear. 
The details depend on whether the ions are fixed or mobile; in either case the 
physical picture corresponds to the single wave model originally proposed by 
0"Neil, Winfrey, and Malmberg for the interaction of a cold weak beam with 
a plasma of fixed ions. 
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I. INTRODUCTION 



Recently, we have studied the collisionless nonlinear evolution of an unstable mode; first 
in a single component Vlasov plasma with fixed ions and then more generally in a multi- 
species Vlasov plasma. [[J- The asymptotic features of the problem in the limit of weak 
instability, i.e. 7 — > + where 7 is the linear growth rate, were our principal focus. The 
main tool has been expansions for the amplitude equation and the distribution functions; 
in particular the asymptotic structure of these expansions. Coefficients of both expansions 
develop singularities as 7 — > + and these singularities reveal the asymptotic features of the 
amplitude equation, distribution functions and electric field. 

The amplitude equation describes the evolution on the unstable manifold of the equi- 
librium and a key conclusion of our previous paper established the scaling behavior of this 
system. || By setting A{t) = ^r{jt) exp(— iO(t)), the resulting equations for r(r) (r = yt) 
and 8(t) were free of singular behavior as 7 — > + , provided the exponent (3 was suitably 
chosen. The correct choice turned out to depend on the model under consideration: if ion 
masses were finite, then typically (3 = 5/2 unless the ion distributions happened to be flat 
at the phase velocity of the linear wave. In the limit of fixed ions (mi — > 00) or when the 
ion distributions are flat at the resonant velocity, the exponent drops to (3 = 1/2. 

In this paper we apply these results for (3 to control and interpret the singularities that 
arise in the expansions of the distribution functions. This study illuminates in detail the 
asymptotic structure of both the distributions and the electric field. In particular, we find 
that the electric field is essentially monochromatic at the wavelength of the linear mode with 
a nonlinear time dependence. Outside the resonant region, the distributions are described 
by the linear eigenfunction and in the resonant region they have a nonlinear structure. The 
details depend on whether the ions are fixed or mobile, but in either case this physical picture 
is well known from the "single wave model" proposed by O'Neil, Winfrey and Malmberg for 
the interaction of cold weak beam with a neutral plasma of mobile electrons and infinitely 
massive ions.. [|J Their work supplied a model of the self-consistent Vlasov problem that 
has proven useful to many researchers in the subsequent years. 0- |12|] Our conclusions 
generalize this useful simplified picture to a general electrostatic instability arising in an 
unmagnetized multi-species Vlasov plasma. As this paper was being completed, we learned 
of the interesting recent work by del-Castillo-Negrete who has given an different derivation 
of the single wave picture using matched asymptotic methods to treat the resonant and 



non-resonant particles. JT3| As in the original work of O'Neil et al, del-Castillo-Negrete 



allows only mobile electrons and moreover restricts attention to instabilities associated with 



so-called "inflection point modes". |l^Jl5| 



In the remainder of this introduction we review our notation and in section [TJ we sum- 
marize the needed conclusions of Ref regarding the singularities of the expansions. The 
third section applies these conclusions to the distributions and electric field, and section [TV] 
contains a final discussion. 



A. Notation 

Our notation follows Ref || ; we consider a one- dimensional, multi-species Vlasov plasma 
defined by 
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(1) 



dE 

dx 



/oo 
dvF {s) (x,v,t). 
-oo 



(2) 



Here x, t and v are measured in units of u/u e , and u, respectively, where a is a 
chosen velocity scale and uj 2 = 4ne 2 n e /m e . The plasma length is L with periodic boundary 
conditions and we adopt the normalization 



L/2 
-L/2 



dx 



dvF {s \x,v,t) 



z s n t 



n. 



(3) 



where q s = ez s is the charge of species s and k^ s ' = q s m e /em s . Note that = —1 for 
electrons and that the normalization (^) for negative species makes the distribution function 
negative. 

Let Fq(v) and f(x,v,t) denote the multi-component fields for the equilibrium and per- 
turbation respectively and k the matrix of mass ratios, 



// (S1) \ 



( Ft ] \ 



V 



K = 



( - ..^ 

k (s2) ■•■ 



/ 



V 



V : J 

then the system ([!]) - @ can be concisely expressed as 
where the linear operator is defined by 

oo 

(Lifi)(v) 



(4) 



/ 



l=— oo 



-il 



1 = 

vfi(v) + k ■ V i(v) Z s I-oo dv> ft\v')} I + 0, 



(5) 

(6) 
(7) 



with T)i(y) = —d v F /l 2 , and the nonlinear operator Af is 

OO OO ' / Q f \ 

Ojm-l 



M(f) = E j mx E 



l=— oo 



I 



K 



dv 



E / dv'fi s '\v'). 



(8) 



In the spatial Fourier expansion (^|), I denotes an integer multiple of the basic wavenumber 
27r/L, and the primed summation in (^) omits the I = term. The notation k ■ r)i(v) 
or k ■ d v f m _i denotes matrix multiplication. For two multi-component fields of (x,v), e.g. 
B = (E>( Sl \ B( S2 \ B (S3 \ ..) and D = (D^ Sl \ ..), we define an inner product by 



r L/2 f oo pL/2 

(B,D)=y dx dvB (s) (x,v)*D (s) (x,v) = / dx<B,D> 

„ J -L/2 J-oo J -L/2 



L/2 
-L/2 



(9) 
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where 



<B,D> dvB^\x,v)*D {s) (x,v). (10) 



The spectral theory for C is well established and the facts needed for our analysis are 
easily summarized. The eigenvalues A = —ilz of C are determined by the roots Ai(z) = of 
the "spectral function", 

J-oo V — Z 

If the contour in (|ll|) is replaced by the Landau contour for Im(z) < then we have the 
linear dielectric €i(z); for lm.(z) > 0, Ai(z) and €i(z) are the same function. The eigenvalues 
can be either real or complex depending on the symmetry and shape of the equilibrium. 

Associated with an eigenvalue A = — ilz is the multi-component eigenfunction ty(x,v) = 
e %ix where 

m = (12) 

v — z 

There is also an associated adjoint eigenfunction ^(x, v) = e ip{y)/L satisfying (\I>, ty) = 1 
with 

$i v ) = ~ T77 V*7~ *\- ( 13 ) 

AJ(^)*(t> — z*) 

Note that all components of i/j(v) are the same. The normalization in ( |l3l) assumes that 
the root of Aj(z) is simple and is chosen so that < ip,^ >= 1- The adjoint determines the 
projection of f(x,v,t) onto the eigenvector, and this projection defines the time-dependent 
amplitude of i.e. A(t) = /). 



II. PREVIOUS RESULTS 

The equilibrium F (v) is assumed to support a "single" unstable mode in the sense 
that E u , the unstable subspace for C, is two-dimensional. With translation symmetry and 
periodic boundary conditions, this is the simplest instability problem that can be posed. 
Henceforth, let k denote the wavenumber of this unstable mode that is associated with 
the root Afc(z ) = which we assume to be simple, i.e. A' k (zo) ^ 0. The corresponding 
eigenvector is 

^( x , v ) = e lkx i){v) = e lkx (-J^Hl) . (14) 

V v — ZqJ 

The root z = v p + / k determines the phase velocity v p = u/k and the growth rate 7 of 
the linear mode as the real and imaginary parts of the eigenvalue A = —ikz = 7 — iu. 
Solutions on the unstable manifold have the form 

f u (x,v,t) = \A{t)i){v)e lkx + A*{t)^*{v)e~ ikx ] + H(x,v, A(t), A*(t)) (15) 
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where A(t) = (9, f u ) evolves according to the amplitude equation 

A = \ A + (^,Af(f u )) 

and self-consistency requires H to satisfy 
f)M f)H 

—A + —A* = CH + Af(f u ) - [(*, Af(f u )) * + cc 



subject to the geometric constraints 



dH 



= H(x,v,0,0) = ^(x,v,0,0) 



dH 

dA* 



x, v, 0, 0). 



(16) 



(17) 



(18) 



The translation symmetry of the model (0) provides important constraints on both the 
amplitude equation and the form of H. || For the amplitude equation (|T6"D , the right hand 
side must have the form 



\A + (y,Af(f u )) = Ap(a) 



(19) 



where a = \A\ 2 and p(cx) is an unknown function to be determined from the model. Similarly, 
translational symmetry requires the spatial Fourier components of H to have a special form 



H (v,A,A*) = ah (v,a) 
H k (v, A, A*) = Aa h ± (v, a) 
H mk (v,A,A*) = A m h m (v,a) 



(20) 



for m > 2 



where = Hi*. These results focus our analysis on a set of functions, {p(a),h m (v,a)}, 
which must be determined from the Vlasov equation. 



A. Expansions and singularities 

We study p(a) and {h m {v, a)} via the expansions 



3=1 



i=i 



(21) 



The coefficients pj and h m j are determined by inserting the expansions into ( |17D and ( |T9"D 
and solving at each order of a. The resulting recursion relations are given in Ref 0] and are 
not required for the present discussion. 

The key point is that for both the amplitude equation and the distribution function the 
expansion coefficients develop singularities in the limit 7 — > + . This can be seen explicitly 
by reviewing the calculation of the cubic coefficient p\. From Ref [|J, p\ depends on h and 

^■2,0; 



Pi 



< d v ip, k ■ (h 0t o - h 2 ,o) > + 1 rr < dvi>, « • V > 



(22) 



where r 2j o = / dv /i 2 ,o- The recursion relations determine ho,o and /i 2)0 , 
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m«) = -^i- 77. k ^ (23) 




M«) = aP (^1+^(^)1^), (24) 

and one notes that for 7 > these are smooth functions but there are complex poles at 
zq and Zq that approach the real axis at v — v p in the limit 7 — > + . For ft,2,o all poles lie 
above the real axis, but ft,o,o contains poles above and below the axis and this forces the 
integral < d v ip, k ■ hop > in p\ to diverge as 7 — ► + because of a pinching singularity. For 
similar reasons, the integral < d v ip, k ■ ip* > also diverges but the remaining integrals in pi 
are nonsingular. 

A detailed evaluation of this asymptotic structure in p\ yields the form 



Pi 



' ^(7) -7^(7) + 0(7*) (25) 



, 1 



7 

where C\ and c?i are nonsingular functions of 7 defined by 

*<7) - -iX^)?^"* 1 - «" ,2 > (£ (*) 

where the primed species sum omits the electrons. At 7 = 0, c% has the limit 

ci(0) = --^yE^ (s) (l - - (S)2 )^K) (28) 

which is typically non-zero yielding a 7 -4 singularity for p\. There are at least two special 
cases of interest for which ci(0) = 0; namely, infinitely massive fixed ions = for all 
s^e) and flat ion distributions at the resonant velocity (77^ (v p ) = for all s^e). In such 
cases, the divergence of p\ drops to 7~ 3 . 

Analagous singularities appear also in the higher order coefficients and grow more severe 
although their character remains the same. The higher coefficients h m j exhibit more and 
more poles which approach the linear phase velocity as 7 — > + and these poles generate 
stronger pinching singularities in the higher coefficients pj. An important property of the 
poles in h m j is that they always have the general form (v — a)~ n or (v — a*)~ n with 

a = z + i-f(/k = v p + i-f(( + l)/k (29) 

where ( > is a purely numerical factor, i.e. the poles always lie along the vertical line 
Re(f ) = v p . 

The explicit calculation of higher order coefficients from recursion relations rapidly be- 
comes prohibitively laborious; however, useful bounds on the singularity of the higher order 
coefficients are obtained Ref || using an induction argument. More precisely, we find for 
the amplitude equation 
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lim 7 !/ |pj| < oo 

7^0+ 



(30) 



for j > 1 where v = 5j — 1 in the generic case with Ci(0) 7^ 0, and u = 4j — 1 in the two 
special cases mentioned above, fixed ions or flat ion distributions, with Ci(0) = 0. For the 
coefficients of the distribution function, the induction argument proves, for m > 0, j > 
and m' > 0, 



lim Y m ' j 
7^0+ 



j— 00 ^ ' 



< OO 



(31) 



where /i mj - = J mj - + 1 with J mj - = (2m + 5j — 3) + 4<5 m0 + 55 m) i in the generic case defined 



by ci(0) 7^ 0. For the special cases of fixed ions or flat ion distributions, (|3lD holds with 
exponent fi mJ = J mJ - j - 5 m ,i- 

The first bound determines the scaling exponent for A. When (0), (|T9|), and ( pl|) are 
combined we obtain an amplitude equation, 



A = \A + Y,PM\ 2j A 
3=1 

where each nonlinear term has a singular coefficient and the equation is ill-defined as 7 
The cure is to rescale the amplitude 

A(t) = >yP r(-yt)e- im 



(32) 



(33) 



with j3 = 5/2 for the typical case (ci(0) 7^ 0) and in the special cases with Ci(0) = we 
require (3 = 2. Once this is done, the equations for r(r) and 9(t) are nonsingular in the 
regime of weak growth rates; additional details may be found in Ref 0. 



III. DISTRIBUTION FUNCTIONS AND ELECTRIC FIELD 

The scaling ( p3[ ) of the amplitude has immediate implications for the asymptotic struc- 
ture of the distributions. From fll5|) the Fourier coefficients of / = F — Fq may be written 
in terms of r(r) and 9(t) in (^) and h m 

f o (v,t) = 1 20 r(T) 2 h o (v, 1 ^r 2 ) 

f k (v,t) = 7^r(r)e-^ [^(v) + 7 2/3 r(r) 2 h^v, 7 2 V)] (34) 
/^(^, = 7 2m/3 r(r) m e" ime « /i m (^, 7 2/3 r 2 ) m > 2. 

As 7 — > + , r(r) is an 0(1) quantity, thus the asymptotic features of each Fourier component 
are determined by the explicitly shown factors of 7 and the asymptotic form of the functions 
ip(v) and h m (v , 7 2/ V 2 ). The dependence on h m necessitates a separate consideration of the 
asymptotic behavior for non-resonant and resonant velocities. In the former regime we 
assume the distance from the linear phase velocity satisfies v — v p = 0(1), and the resonant 
regime corresponds to velocities within a neighborhood of v p that scales with the growth 
rate, i.e. v — v p = 0(7). For resonant velocities, the singularities in ip(v) and h m (v , 7 2/ V 2 ) 
come into play and alter the asymptotic features of the distribution function. 
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A. Non-resonant velocities 



For v — v p = 0(1), the functions ip{v) and h m (v, a) are bounded 0(1) quantities (we use 
a = -y 2 / 3 r 2 to emphasize this), and the Fourier components (33) combine to yield 



(F{x,v,t)-F (v))/^ 



r(r)e- im ^(x, v) + cc] + 7^ r(r) 2 h (v, a) 

imkx 



(35) 



+7 2/3 r(r) 3 



+ E [ 7 (2m - 1)/3 e- imeW hm(v, a)e imkx + cc 



m=2 



r(r)e~ im ^(x, v) + ccl + 



(36) 



In words, the non- resonant correction to Fq scales overall as 7^; the leading piece of this 
correction simply has the form of the linear wave *&(x, v) with nonlinear time dependence 
determined by the mode amplitude r(r) exp(— i9{t)). 



B. Resonant velocities 

For v — v p = 0(7), the functions ip(v) and h m (v,a) typically develop singularities as 
7 — > + and these divergences compete with the explicit factors of 7 in (B3J) to determine 
the asymptotic form of the distribution. The analysis is simplified by the fact that all relevant 
singularities are poles of the form described in (|29|), and these may be rewritten as a singular 
factor multiplying a nonsingular function of the rescaled velocity variable u = (v — u p )/7, 
e.g. 

1 = 1 1 = 1 1 (3 7) 

(v-a) n 7 n {{v -VjO/7 -i(C + l)/^) n 7 n ((u - z(C + 1)/A;) n ' 

Once this is done, the functions ^(u) and h m (v,a), expressed in terms of u, may be sub- 
stituted into Q3"3|); the variable w provides a uniform velocity coordinate for the resonant 
region. 

The puzzle is to deduce the correct overall factor of 1/7™ for each function. For h m (v, a) 
we have the integral bound ([H]) which may be rewritten in terms of u 



lim 

7-S-0+ 



-00 v ' 



< 00. (38) 



Since all singularities are poles we know the integrand does not have an integrable singularity, 
hence we conclude that 7 1+ ^ m ^ h m j(v p + ju) defines a nonsingular function of u: 

h m A v P + 7«) = T {l+ ^ y h m j(ui 7)- (39) 

The nonsingular character of h m} j(u,j) can be checked directly for the specific examples in 
(E3I) and (PH), and also verified, in general, from the recursion relations. From the expansion 



of h m (v, a), we thus find 
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h m (v p + 7 «, a) = £ T 2 ^~ (1+Mm ' j) A m j(w, 7)r 2i - (40) 

j=0 

In the generic case with /3 = 5/2 and /x mj - = J mj - + 1 this gives 

■y oo 

h m (v P + 7«, a) = — ^ 7)r 2j (41) 

7 m j= o 

where 

3 m = 

6 m = 1 (ci(0) ^ 0). (42) 

2m — 1 m > 2 

In the special cases with fixed ions or flat distributions, then (3 = 2 and /i mj - = J m ,j—j — S mi x, 
and (fill) holds with exponent 

2 m = 

5 m = <( 4 m = 1 (ci(0) = 0). (43) 

2m - 2 m > 2 

In all cases, we define the nonsingular function h m (u, r 2 , 7) = Y^jLo h m j(u, 7)r 2j and rewrite 

m 



h m (v p + iu,a) = hm{u f> 1 \ (44) 

It is simpler to obtain the corresponding factorization of the eigenfunction; from the 
definition (|14]) we have 

^{v p + 7 M = , 45 

7 \ u — i/k J 

and the only subtlety concerns r)k(v p + ju) which is 0{1) in the generic case (ci(0) 7^ 0) 
and 0(7) in the two special cases with ci(0) = 0. Thus we define the nonsingular function 

i){u,i) by 

^ v+lu) = tkhA ( Cl (0)^0), (46) 

7 

in the generic case, but in the special cases the eigenfunction is itself nonsingular and we 
have 

iP(v p + ju) = </>(«, 7) (ci(0) = 0). (47) 

We are now able to describe the asymptotic structure of the distributions. 
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1. Generic instability: ci(0) ^ 
For the generic case, inserting fl44 ) and ([16|) into (|34"D yields 



[F(x, v p + ju, t) - F (v p + 7«)]/7 3/2 = {r(r) 



-i0(t) ikx 



$(u, 7) + r(r) 2 /ii(w,r 2 ,7) + ccj 



^(m-2)/2 gimfex r ^ r ^m g - 



eo 



m=2 



(4* 



neglecting the subdominant terms this gives 
[F(x, v p + 714, t) - F (t; p + 7 u)] / 7 3/2 = {r(r) e «» 7 ) + r ( r )2 ^{u, r 2 , 7)] + cc} 



(49) 



The generic resonant correction to F , expressed in the velocity coordinate u, scales overall as 
7 3//2 ; the leading term in this correction has the wavelength of the linear wave but the velocity 
dependence, ip(u, 7) +r(r) 2 hi(u, r 2 , 7), is not simply given by the linear eigenfunction. The 
time dependence is determined by the mode amplitude r(r) exp(— i9(t)) but the dependence 
on r is rather complicated. 



2. Special cases: c\ (0) = 

For the special cases, defined by fixed ions or flat ion distributions, we apply (K7|) and 
to d34T) and obtain 



[F(x, v p + 7u, t) - F (v p + 7 m)]/7 2 = {r(r) 



-id(t) ikx 



7) + r(r) 2 hi(u,r 2 ,^y) + cc} 



<r) 2 ^o( U , r 2 , 7) + £ [e mte r (r) m e - imfl «> ^(u, r 2 , 7) + cc 



771=2 



(50) 



This is a qualitatively different structure in contrast to (HS); now the resonant correction 
is £?(7 2 ) and a// wavelengths are present at leading order. Thus the spatial dependence is 
very rich and bears no special relation to the linear instability; a similar observation holds 
for the dependence on velocity. 



C. Electric field 



The Fourier components of E are given by fl34|) and Poisson's equation 



ikE k (t) = 7 /3 r(r)e 



-i0(t) 



J —00 



(51) 



/OO 
*>X>m(*>,<r) m>2. 
-00 a 

Bounds on the asymptotic form of the integrals can be inferred from the expansion h r 



J2j h m j(j 3 and the bound (|3lD on the integrals of h m j (for m! — 0). The details depend on 



whether we consider the generic instability or the special cases. 
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1. Generic instability: ci(0) ^ 



From (|3~ID we find 



lim 7 a " 

7^0+ 



dvJ2h£H 



< oo 



with 



I 2 


m 


= 






m 


= 1 


(ci(O)^O) 


[ 2m 


— 2m 


> 2 





Hence, with /3 = 5/2, the generic components are 
ikE k (t) =7 5 / 2 r(r)e- ifl(t) 



imkE mk {t) = 7 3m+2 r(r) m e _imfl W ^-2 



l + r(r) 2 7 5 /°° tfo;X;M a) 



The asymptotic electric field is 



£7(a;,t) _ 1 
" Jfe 



7 



5/2 



-ir(r)e 



-i6{t) 



l + r(r) 2 7 5 /°° cfo^/4 s) (w,a) 



e ikx + cc 



(52) 



(53) 



(54) 



m > 2. 



0(7 11/2 ); (55) 



clearly is dominated by the wavenumber of the unstable mode with an overall scaling of 
7 5 / 2 . The term 7 s / dv J2 S ^1 is treated as an (9(1) contribution in light of the estimate ( |52"D 
above. 



2. Special cases: c\ (0) = 



For instabilities with fixed ions or flat ion distributions, we have (3 = 2 and fj, m j 
1 in (|31|) ; applying this bound to the series h m = Ej h m ja j yields 



lim 7 a " 

7^0+ 



< 00 



with 



a, 



1 


m 


= 




.3 


m 


= 1 


( Cl (0) = 0) 


2m 


— 3m 


> 2 





Now the general expressions for the components reduce to 



ikE k {t) = 7 2 r(r)e 



-ie(t) 



l + 7r(r) 2 7 3 / dv^h^i 

J —OO 



/OO 
-00 „ 



(56) 



(57) 



(58) 



m > 2, 
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and the asymptotic electric field has the form 



E(x,t) 
7 2 k 



i {-ir(r)e- ie ^ [1 + 0( 7 )] e ifcx + cc} + C( 7 5 ). (59) 



The overall scaling is now the well known 7 2 or "trapping scaling" and the leading term has 
a much simpler structure. Again we find the wavenumber k of the linear instability; however 
now the time dependence is simply given by the mode amplitude r(r) exp(— i9(t)). 

IV. DISCUSSION 

The single wave model, derived originally by O'Neil, Winfrey and Malmberg, described 
the interaction of a cold electron beam interacting with a plasma of mobile electrons and fixed 
ions. In their problem, the infinite extent of the plasma allowed for continuous wavenumbers 
and the dispersion relation for a cold beam was required to select a single wavenumber 
corresponding to the maximum growth rate. This wavenumber characterizes the electric 
field whose nonlinear time development results from the coupling to resonant particles. The 
nonresonant plasma simply provides a linear dielectric which supports the wave. 

By contrast, we pose a more general problem, allowing for multiple mobile species and not 
restricting the type of electrostatic instability, but for a finite plasma with periodic boundary 
conditions. Within this setting, we consider equilibria supporting a single unstable mode 
and derive the resulting equations for the electric field and distributions in the limit of 
weak instability. In this asymptotic limit, the physical picture of the original single wave 
model emerges quite generally. The monochromatic electric field is coupled to the resonant 
particles and evolves nonlinearly while the nonresonant particles show only a linear response 
to the electric field. 

The amplitude expansions, whose singularity structure form the basis of our analysis, 
do not provide a practical tool for solving the single wave model. For this purpose, it is 
more convenient to assume the simplifications of the single wave picture and derive model 
equations directly from the original Vlasov theory. This development will be presented in a 



forthcoming paper. [|XC 
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